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ABSTRACT 

The velocity distribution in the laminar boundary layers at the interface of 
two liquids in relative rectilinear motion, the thickness of the layers, and the 
stress at the interface are determined. Numerical results are given for nine cases 
of liquids in contact, including identical liquids and liquids with varying degrees 
of dissimilarity in characteristics. The evaluation of the desired quantities is 
based on PrandtPs boundary-layer theory, and is carried out by a method of suc- 
cessive approximations. The numerical results are those given by the second 
approximation. 
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I. INTRODUCTION 

When a liquid of given density is flowing over a liquid of greater 
density, the latter being at rest, three successive forms of the inter- 
face are discernible as the velocity of the moving liquid is increased. 
To take a concrete example, let us assume that the flow takes place 
in a closed rectangular channel with a lower pool; the upper liquid 
is fresh water and the lower liquid is a solution of salt or sugar in 
water. First, at small velocities, the surface of separation is sharp 
and distinct, indicating, since the indices of refraction of the liquids 
are different, that the densities are discontinuous at the interface 
and the flow at the interface is laminar. The flow in the central 
part of the upper liquid, on the other hand, may be laminar or turbu- 
lent, depending on the size of the channel. Secondly, if the velocity 
of the upper liquid is gradually increased, at some critical velocity 
the smoothness of the separating surface first disappears, and then the 
surface becomes covered with surface waves traveling in the direction 
of the current and with a velocity slightly smaller than the velocity 
of the upper current. When the waves first appear, they are long- 
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crested, that is, the length of the crests normal to the direction of 
motion is greater than the spacing between the crests. The crest 
lines are parabolic, concave upstream, and highest at the center of 
the channel, but flattening out toward the channel wall. At this 
stage of the flow — and this perhaps is significant — the waves are 
stable and show no tendency to break. The waves travel with 
practically no deformation, except the tendency to grow slightly 
larger during their motion. Finally, when the velocity of the upper 
current exceeds the critical velocity, the waves become sharp-crested 
and the crests shorter. The waves no longer are stable; that is, 
portions of the crests break away from the waves and are thrown into 
the upper current. These portions move forward and upward. In 
fact, this is the manner by which the mixing of the lower liquid with 
the moving liquid is brought about. When the velocity of the 
current is increased, still further within certain limits, the rate of 
mixing increases. At this stage of the flow, and this also is perhaps 
significant, the wave length is not affected by changes of velocity. 
These are the main results of observations made in a laboratory channel 
-here the depth of the upper current is not over 10 centimeters and 
the length of the channel is not over 10 meters. 

A detailed explanation of these phenomena is lacking. A theoretical 
approach supplementing experimental knowledge should prove help- 
ful. A theoretical investigation should begin with a study of the 
stability of the interface on the basis of the actual conditions in the 
neighborhood of the interface. One of these conditions is a discon- 
tinuity of density at the surface of separation. Another is a pair of 
laminar boundary layers on the two sides of the surface of separation. 
In this paper the thicknesses of these boundary layers and the velocity 
distributions within them are determined by means of successive 
approximations, on the assumption that the interface is a smooth, 
plane surface. 

II. FORMULATION OF THE PROBLEM 

A liquid of indefinite height and having initially a uniform velocity, 
U, meets and flows over a still liquid of indefinite depth. In general, 
the two liquids have different physical characteristics; that is, their 
viscosities and densities are not the same. It will be assumed that 
the lower liquid has the greater density. If the pressure is everywhere 
hydrostatic, the interface or surface of separation is a horizontal 
plane. When U is sufficiently small the narrow region on both sides 
of the interface is a region of viscous laminar flow. Thus, the acti- 
vating forces at the point where the moving liquid comes first in con- 
tact with the lower liquid are tangential and act in the horizontal 
direction; there is no normal force acting on the interface and hence 
the interface remains horizontal. Every particle of liquid on the 
interface moves with a constant velocity, u . The thickness of the 
viscous boundary layers in the two liquids is initially zero and increases 
with the length downstream (see fig. 1). In the upper layer the 
velocity in a vertical section increases from u to U; in the lower 
layer it decreases from u to 0. It is proposed to evaluate the variation 
of the thickness of the laminar layer with the length of the interface, 
the variation of velocity in the two layers, the shearing stress at the 
interface and the dependence of these quantities on the physical 
characteristics of the two liquids. 
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Figure 1 . — Laminar boundary layers at the interface of two liquids. 
III. PRANDTL'S BOUNDARY-LAYER EQUATION 

In order to avoid negative values of ordinates two distinct sets of 
axes will be used, x, y for the upper liquid and x' ', y f for the lower 
liquid. See figure 1. The a:-axes of the two sets coincide and lie 
in the plane of separation of the two liquids, the positive sense being 
in the direction of L . The point at which the two liquids first come 
into contact will be taken as the origin. In the upper liquid, positive 
y is directed upward; in the lower liquid, positive y f is directed down- 
ward. The density, the viscosity, and the kinematic viscosity of the 
upper liquid will be denoted by p, /x and v, respectively, and the 
corresponding properties of the lower liquid by p', m', and /, respec- 
tively. The parameter, r, defined as 



^=V[/p7mp], 



(i) 



is of special significance in this analysis, as will be seen subsequently. 
The basis of these computations is the boundary-layer equation 
of Prandtl for the two-dimensional flow of an incompressible fluid [l]. 1 
This equation with the law of hydrostatic pressure and the equation 
of continuity, when expressed in terms of the upper liquid, becomes 






dx 



d 2 u 



I_dp 
p dx' 



1 dp. 



0, 



du dv 
dx^dy ' 



(2) 
(3) 
(4) 



where u and v are the components of velocity in the x and y direc- 
tions, respectively, and p is the hydrostatic pressure. 

» The figures in brackets indicate the literature references at the end of this paper. 



306 Journal of Research of the National Bureau of Standards 

In view of the two-dimensional equation of continuity of an incom- 
pressible fluid, eq 4, the stream function, $, may be introduced: 

U —W V =Zx' (5) 

Since, in the present case dp/dx=0, eq 2, in terms of the stream 
function, becomes 

by dxby bx dy 2 ~~ v by 3 ^ ' 

This can be changed into an ordinary differential equation by choos- 
ing appropriate new variables. The physical basis for the ordinary 
differential equation is the similarity of flow in two vertical sections. 
A necessary condition for the similarity of flow is that the velocity 
at the interface u and the velocity at large values of y be constant. 
These conditions are fulfilled if the pressures are hydrostatic, which 
has already been assumed. Adopting a characteristic length 5, 
which is to be regarded as a function of x only, the new dimensionless 
variables for y and the stream function are 

v=ny/t, (7) 

and 

H( v ) = -+/U8, (8) 

where n is a dimensionless numerical constant as yet not specified. 
These variables are dimensionless. In terms of these, eq 6 now 
becomes 

TJ d8 „ d 2 H . v d?H h m . 

u Tx H W +n -sW= - (9) 

Thus, if we select 5 to satisfy the relation 

dS v 



or 



h dx =n V 



S 2 =2n % (10) 



then, from eq 9, we have 



U' 



3?+*£-* <") 



which is Blasius' equation for the laminar boundary layer when the 
pressure is independent of x [2]. 

In eq 10 the dependence of 8 on x is given except for the numerical 
value of n. We determine this by considering the velocity at large 
values of y. Now, 

dy drj 

hence if we specify that dH/dr]=l for large values of 77, and this is a 
boundary condition for the upper fluid, we must then take n=\. 
These values we now adopt. 



" ' " ■ , ■ ' 



Flow at Liquid Interfaces 307 



Summarizing, we have for the upper liquid 



(P7J d 2 TT 
i«=ff+/4f=0, (12) 



where 



H=— jj-jrj^-, (13) 



and 



=Vt- m> 



A similar set of equations may be established for the lower liquid by 
using primes to distinguish the variables, 

,7377/ ^277/ 

W)^+3<^-0, (15) 



where 



H '=-m»i'=w' ^ 



and 



'=V 2 # (,7 > 



for the upper liquid; and 



The actual velocities and their first derivatives in terms of the 
dimensionless variables are 

«=tff, (18) 

bu_Ud 2 H . . 

dy~ 8 d v 2 ' Uy; 

u>=U<%jr, (20) 

W = Ud 2 H> 

dy' 8' drj 
for the lower liquid. 

For the simultaneous solution of the two basic differential equa- 
tions, eq 12 and 15, six boundary conditions are needed. These we 
now consider. In the upper liquid, u equals U for y=&>', in the lower 
liquid, u' vanishes for y' = <x> . Accordingly, 

ff=i, ,=»; (22) 



7 -Jir-* (2D 



and 



fr=<w=». (23) 
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The stream function, \[/ y may be interpreted as the flux per unit width 
between a point y of the upper liquid and the interface; and $' as the 
flux per unit width between a point y' of the lower liquid and the 
interface. As we have supposed that the interface is the horizontal 
plane y=y' = 0, we must have 

ff=0, 7?-0; (24) 

and 

H'= o, tj' = 0. (25) 

The velocities are continuous at the interface; hence 

drj drj 
Again, the shearing forces are continuous at the interface ; hence 



7> r)=v' = 0- (26) 



du ,du' , 

m d 2 H / dm' 

or, using eq 1, 

d*H d*H' , n ,„. 

These are the six boundary conditions of the problem. 

IV. APPROXIMATE EXPRESSIONS FOR THE STREAM 

FUNCTIONS 

To solve the differential equations, eq 12 and 15, we shall resort to 
a method of approximations which is a modification of the method 
used by Pohlhausen for the solution of the Blasius plate problem [3]. 

In this method certain approximate boundary conditions for large 
values of 77 and 77' are required. First, consider the upper moving 
liquid. Select the value 77 s of 77 such that 

l-^£;XnZn,, (28) 

where e is a small quantity, say 1/100 or any other smaller fraction. 
Denote the value of H corresponding to rj s by H s . Now, writing 
eq 12 in the form 

drf ~ H d V 3 ' ( " V) 

and integrating between 7? = ij s and »j=°°, we have 



\dvh~ J, 8 Hd v z dr >- 



Integrating the right-hand member by parts, and neglecting the 
term multiplied by e, and making an obvious substitution from eq 29, 
using eq 28, we obtain 



(dH\ _ Wd 2 H\ f°° 1 d?H , 
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Treating the integral in the right-hand member again in the same 
manner, we obtain 

,J d M\ -(*K\ (1 L\ * n^, 

UvK WA.W. II}) 6 J Vs H>d/ dv - 

If € is sufficiently small, that is, if H s is sufficiently large, the integral 
in the latter equation may be discarded, leaving 



\dvA, WAU* H?) 



(30) 



This is the approximate relation which exists between the first and 
second derivatives of H when 77=77*. 

If we now suppose that H s is so large that 1 can be neglected com- 
pared to the square of H s , then, in view of eq 28, 29, and 30, we must 
have 

dH \ 

df\ 

dm 
am 

drj z 



= 1-6, 
= €#„ 






-jjr^-tHs 2 , y=y s . 



(31) 



Consider now the lower still liquid. We select the value 77/ of 77', 
such that 



dJV 

f £«', if •**', 



(32) 



where e has the same general meaning as e. Denote the value of 
H' corresponding to 77/ by H/. For values of 77' larger than 77/, eq 
15 may be written in the form 



dm' ,d 2 H' 

dr,'* ~ Hs drj' 2 ' 



Obviouslv, then 



z *n* 



dH' , 

dr>' ~ € > 

jri — — e tl t , 77 —77s , 

j /3 —etls, 77 —77s . 



(33) 



Now, in the method employed here, the end conditions for the 
large values of 77 and 77 ' will be taken from eq 3 1 and 33 after assigning 
definite values to e and e' '. The values assigned to e and e' have a 
bearing on the accuracy of the approximation and also govern the 
course of the computations. If e and e' are assigned values different 
from zero, the accuracy of the approximation is increased, but the anal- 
ysis becomes more complicated. In order to simplify computations, we 
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shall take e=e' = and adopt the following end conditions. For the 
upper liquid, we have 

r^W^'W^ ' 11 ^' (34) 

and for the lower liquid, 

dH' n d?H> cPH' , , 

W^ ' ~d7*= 0> W =0 > " = " s • (3i>) 

These now replace the conditions implied by eq 22 and 23. With 
this choice of the end conditions, rj s and r// have simple geometric 
interpretations. Let y s and y s ' be the values of the distance from the 
interface corresponding to rj s and 17/, respectively; that is, let 



y,=r,,j(2vx/U), (36) 

and 

y/=n/V(2^/C7), (37) 

then y s and y/ are the thicknesses of the two laminar boundary 
layers at the interface. 

We shall also need the following end conditions from eq 24, 25, 26, 
for the upper liquid, 

dtt 

ff=0, ~—a , 77=0, (38) 

and for the lower liquid, 

#'=0, 57=< r>'=0, (39) 

where 

a =a '=u IU. 

We are now ready to develop the general expressions for the stream 
functions H and B.' . Since from eq 15 and 34, d 3 H/dt] d vanishes at 
rj=0 and rj=r} s , then it must be possible to express this function of H 
as a summation of polynomials where each polynomial vanishes for 
r;=0 and 7j=rj s . Suitable polynomials may be selected in various 
ways. In order to extend Pohlhausen's computations [3], we select 
the following oscillating algebraic functions as the expansion poly- 
nomials for this analysis: 

Ni= — a(a— 1) = a— a 2 , 

N 2 =-+a(a-l)(2a-l)=2a 3 -3a 2 +a, 

N 3= - a ( a -l)(Sa-l)(3a-2)=2a-lla 2 +18a~9a\ 

and, in general, 

N n ^(-l) n a(oi-l)(na-l)(jia~2)(na-3)( . . . )(na-n+l), 

where 

a=ri/r) s . 
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Using only the first three functions, we put 

^■=^(-6^+32*^2-30*8^8). (40) 

Here a 1} rj s , k 2 , and k 2 are constants to be determined. Integrating 
once with respect to rj, we have 

d 2 H 

-j-2-=ai(J l +k 2 j2+h^3) 9 (41) 

where J i} J 2 and J z are functions of a: 

J 1 ^l_3a 2 +2o: 3 , 
J 2 =16a 2 -S2a 3 +Ua\ 
J 3 =l-30a 2 + 110a 3 - 135a 4 + 54a 5 . 

It must be remarked that J x is derived from iVi by integrating with 
respect to r\ and selecting the constant of integration so that Ji = 
for a=l, and similarly for J 2 and J 3 . Equation 41 thus obviously 
satisfies the second end conditions in eq 34. Integrating eq 41 with 
respect to rj, we have 

^=<h+aMG 1 +hG 2 +hG 3 ), (42) 

where G lt G 2 , and G 3 are functions of a: 

G^a-oP+^a*, 

<? 2 =^a 3 -8a 4 +^ct 5 , 

£ 3 = a-10a 3 +ffa 4 -27a 5 +9a 6 . 

In this case, G x is derived from J x hy integrating the latter with re- 
spect to 17 and selecting the constant of integration so that (ti=0 
when j?=0, and similarly for G 2 and G 3 . The constant a is intro- 
duced in order to satisfy the second end conditions in eq 38. Inte- 
grating eq 42 with respect to 77, we have 

11= aorj^+a^ (F t + k 2 F 2 +k 3 F 3 ) , (43) 

where F , F 1; F 2 , and F 3 are functions of a: 

F Q =a, 



4 
: 3 a 5" ' 15* 



F 2 =^-fa 6 +^y 



F 3 =|a 2 -|a 4 4-^a 5 -|a 6 4-|a^ 
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where F x is derived from G x by integrating the latter with respect to 
7j and selecting the constant of integration so that 2^ = when 77=0, 
and similarly for F , F 2 , and F 3 . It is obvious that the first of the 
end conditions in eq 38 is satisfied. The functions N, J, G, and F 
are tabulated in tables 1, 2, 3, and 4. 



Table 1. — Values of the N functions 


a 


Ni(a) 


iV 2 (a)X10 


A T 3(a)Xl0 


0.00 


0. 0000 


0. 0000 


0. 0000 


.05 


.0475 


.4275 


.7030 


.10 


.0900 


.7200 


1.0710 


.15 


.1275 


.0925 


1. 0869 


.20 


.1600 


.9600 


0. 8960 


.25 


.1875 


.9375 


.5859 


.30 


.2100 


.8400 


.2310 


.35 


.2275 


. 6825 


-. 1081 


.40 


.2400 


.4800 


-.3840 


.45 


.2475 


.2475 


-. 5631 


.50 


.2500 


-.0000 


-.6250 


.55 


.2475 


-.2475 


-.5631 


.60 


.2400 


-.4800 


-. 3840 


.65 


.2275 


-.6825 


-. 1081 


.70 


.2100 


-.8400 


.2310 


.75 


. 1875 


-. 9375 


.5859 


.80 


.1600 


-.9600 


.8960 


.85 


.1275 


-. 8925 


1. 0869 


.90 


. 0900 


-.7200 


1. 0710 


.95 


.0475 


-.4275 


0. 7030 


1.00 


.0000 


-.0000 


.0000 



Table 2. — Values of the J functions 



a 


</i(a) 


Ji{a) 


./3(a) 


0.00 


1.0000 


0. 0000 


1. 0000 


.05 


0. 9929 


.0344 


0. 9451 


.10 


.9720 


.1296 


.7971 


.15 


.9393 


.2632 


.6312 


.20 


.8960 


.4096 


.4842 


.25 


.8437 


.5624 


.3690 


.30 


.7840 


.7056 


.3075 


.35 


.7183 


.8280 


.2979 


.40 


.6480 


.9216 


.3366 


.45 


.5747 


.9800 


.4105 


.50 


.5000 


1. 0000 


.5001 


.55 


.4253 


0. 9800 


.5904 


.60 


.3520 


.9216 


.6633 


.65 


.2817 


.8280 


.7008 


.70 


.2160 


.7056 


.6942 


.75 


.1563 


.5624 


.6309 


.80 


.1040 


.4096 


.5280 


.85 


.0607 


.2632 


.3675 


.90 


.0280 


.1296 


.2037 


.95 


.0073 


. 0344 


.0615 


1.00 


.0000 


.0000 


.0000 



p. 
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a 


Gi(«) 


0i (a) 


Gs(«) 


0.00 


0. 0000 


0. 0000 


0. 0000 


.05 


.0499 


.0011 


. 0480 


.10 


.0990 


.0045 


.0927 


.15 


.1469 


.0144 


.1270 


.20 


.1928 


.0299 


.1569 


.25 


.2363 


.0552 


.1760 


.30 


.2770 


.0869 


.1942 


.35 


.3146 


.1254 


.2093 


.40 


.3488 


.1691 


.2255 


.45 


.3793 


.2168 


.2444 


.50 


.4062 


.2669 


.2641 


.55 


.4294 


.3165 


.2932 


.60 


.4488 


. 3640 


.3244 


.65 


.4646 


.4077 


.3590 


.70 


.4771 


.4464 


.3937 


.75 


. 4863 


.4781 


.4272 


.80 


. 4928 


.5024 


.4559 


.85 


. 4969 


.5191 


.4783 


.90 


. 4990 


.5288 


.4925 


.95 


. 4999 


. 5326 


. 4999 


1.00 


.5000 


.5333 


.5000 



Table 


4. — Values of the F functions 


a 


Fi(a) 


Fi (a) 


Ft (a) 


0. 00 


0. 0000 


0.0000 


0. 0000 


.05 


.0012 


.0000 


.0012 


.10 


. 0050 


.0001 . 0047 


.15 


.0112 


. 0006 


.0106 


.20 


.0196 


. 0016 


.0176 


.25 


. 0303 


.0037 


. 0261 


.30 


.0432 


.0073 


. 0351 


.35 


. 0580 


.0127 


.0449 


.40 


.0746 


. 0200 


. 0556 


.45 


.0928 


.0297 


.0692 


.50 


.1125 


.0415 


. 0807 


.55 


.1334 


. 0562 


.0995 


.60 


. 1554 


.0733 


.1096 


.65 


.1782 


. 0925 


.1267 


.70 


. 2018 


.1138 


.1457 


.75 


.2258 


.1370 


.1661 


.80 


.2505 


.1618 


.1881 


.85 


.2753 


.1874 


.2118 


.90 


.3001 


.2134 


.2357 


.95 


.3750 


.2398 


.2608 


1.00 


.3500 


.2666 


.2857 



In these expressions there appear two types of constants. The 
constants k 2 and k 3 are of one type and may be referred to as the ex- 
pansion constants, since these are introduced because of N 2 and N 3 . 
The constants a , a h and tj s are of the other type and may be referred 
to as the basic constants, since a relates to the interfacial velocity 
u , a h to the local shear and v s to the thickness of the laminar boundary 
layer in the upper liquid. 

Reasoning in a similar manner, we have for the lower, still liquid, 



cPH' «! 



-,(-6iV 1 + 32i 2 W 2 -30& 3 W 3 ), 



drj'* 
d 2 H' 

dH' 



and 



dr,' 



=a o '+aiV«7i+fe'0 ! .+Jfca'<7»), 



H'=a,'r l ,'F a +a l ' v P{F 1 +h'F i +h'F z ), 



(44) 



(45) 
(46) 



(47) 
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where N, J, G, and F are the same functions as before but now depend 
on the independent variable of a'=n'lnf. Again, k 2 ' and k/ are the 
expansion constants, and a ', a/, and 77/ are the basic constants for 
the laminar boundary layer in the lower, still liquid. 

Equations 43 and 47 are the approximate expressions for the stream 
functions H and H' . The unknown elements are the basic and expan- 
sion constants, which we now proceed to determine. 

V. EVALUATION OF THE CONSTANTS 

There are a total of six basic constants and four expansion constants 
in the expressions for H and H' , and thus 10 relations are required for 
their determination. Four such relations are obtained immediately 
by considering the first of the end conditions in eq 34 and 35 and the 
continuity conditions, eq 26 and 27. These give, first, 

a +a 1 r} s A =l, (48) 

A =l(l+^k 2 +hJ, (49) 

Oo'+^V^o'^O, (50) 

A'^i+^h'+h')' (5i) 

a =a '; (52) 

a^-ra^M; (53) 

M=(l+h')l(l+h), tf=(/pO/(MP). (54) 

The remaining six relations may be obtained from eq 12 and 15 by 
using one of two methods. 

In the first method the left-hand members of eq 12 and 15 are mul- 
tiplied successively by JVi, N 2 , and Nz, and integrated with respect to 
a or a! between the limits and 1 and set equal to zero, since L is 
identically equal to zero. That is, putting 



where 

second, 
where 

third, 

and fourth, 

where 



and 






—— 



'- " ■' »'■' ,/ » "■ ' " 
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the desired six relations are 

f 1 LN t da=Q, (i=l,2,3); 



315 



and 



(57, 58, 59) 

(60,61,62) 

A simpler method would be to use ft, Q 2 , and Q d instead of N u N 2) 
and N 3 , where 



flL'N t da' = 9 (i=l,2, 3). 



Qi= + 1, Ogaora^l; 




Q 2 = + l, 0<aor a'<l/2; 
Q 2 = — 1, l/2^aor c/jSl; 




8 = + l, 0<aor a'gl/3; 
Q 3 = — 1, l/3<aor a'<2/3; 
Q 3 = + l, 2/3^ a or a'gl, 




The second method, which is the more familiar one, 


utilizes the 


following six relations [4]: 




f 1 o Lda=0; 


(63) 


flL'da=0; 


(64) 


dL 


(65) 


f-». ^ 


(66) 




(67) 



dL' 



drj 



7=0, u' = 0. 



(68) 



The constants of the present problem have been determined, by- 
using both of the above methods, and the results agreed to the same 
order of accuracy. The results obtained from eq 63 to 68 will be 
given here. 

From eq 40, 41, and 43, we find from eq 55, 

L(v)=- (-6M + 32i 2 2V 2 -30&8# 3 ) 

Vs 

+ (whs (JiFo + k 2 J 2 F + faJ 3 F ) 

+a l %\J 1 F 1 +fa[J 2 F 1 +FM+fa[J i F l +F 9 J l ]) 

+a 1 %\k 2 2 J 2 F 2 +hUJ2F- i +F 2 J 3 ]+k^F z ). 
Substituting in eq 63 and integrating, we find that 
a x 2 = A^a^ 2 + A^rjs 3 , 



where 
and 



-Ai— (<ki+<k2k 2 +awkz)/(l+k z ) } 

A 2 = (an + a l2 k 2 + a lz k 3 + a 22 k 2 2 + a 2Z k 2 fa + a Z3 fa 2 ) / (1 + fa) . 



(69) 

(70) 
(71) 



and 
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The values of the coefficients with the double subscripts occurring 
in eq 70 and 71 are 

Ooi = 0.15000, a u = 0.029365, a 22 = 0.0308, 

a 02 = 0.26667, a 12 —0.07619, a 23 =0.0758, 

a 03 =0.21425, a 13 — 0.07625, a 33 = 0.0400. 

By treating eq 64 in a similar manner, we find that 

where 

A/=(a 0l + a O2 k 2 ' + a zh')l(l+h'), (73) 

LU 

^ 2 '=(a n + a 12 * 2 '+a 13 fc 3 ' + ^^ (74) 

Since from eq 55 

dL_dm dH d 2 H am 
diTdif + di dii^ tJ - drf' 

and 

dA\ dN 2 dN 3 9 n n 

It" 1 ' ~da =1 > da= 2 > ^°' a = > 

dNi 1 rf2V 2 - dN 9 ' 1 

— y— = — 1, -j — =1, —j — = —2. 77 = 77^, a=l, 

we find from eq 65 

^ (6+32*2+60*3) = 0, (75) 

and from eq 66 

^(-6 + 32* 2 -60t 3 ) = -a a 1 (l+* 3 ). (76) 

Vs 

Adding eq 75 to eq 76, we have 

64t 2 =-a r; s 2 (l + * 3 ). (77) 

Since i? 8 is finite and a x does not vanish, we have, from eq 75, 

60& 3 =-6-32* 2 . (78) 

We have from eq 77, using eq 48, the relation 

64fc 3 =- 4ffio(1 ffl 7 ao)2 [(l+^)/(l+^ 2 + Ar 3 ) 2 ) (79) 

where t) s does not appear. It may be emphasized that eq 78 and 79 
are two relations expressing the expansion constants k 2 and * 3 in 
terms of the basic constants a and a x . 
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In a similar manner we obtain from eq 67 and 68 the two relations 

60£ 3 '=-6-32£ 2 ', (80) 

6W = -^ 3 [(l+^0/(l + ^+^) 2 } (81) 

Equations 48, 50, 52, 53, 69, 72, 78, 79, 80, and 81 are the desired 10 
relations for determining the 10 constants appearing in the expressions 
for H and H f given by eq 43 and 47, respectively. In principle the 
determination of the constants in a direct manner is possible, but 
since this is tedious, it may be carried out by successive approxima- 
tions. The successive approximations to the values of the constants 
give the successive approximations to the present boundary-layer 
problem. 

VI. FIRST APPROXIMATE SOLUTION FOR THE LAMINAR 
BOUNDARY LAYERS 

In the first approximate solution it is assumed that the expansion 
constants k 2 , & 3 , k 2 , and k z f are zero. 
Consider the relations 

Oo+aii?*Ao=l, (48a) 

ao' + aiVA'^O, (50a) 

ao=ao', (52a) 

a^-ro/A/, [r=(MV)/(MP)], (53a) 

a 1 2 =A l a a 1 2 rj s 2 +A 2 a 1 :i r) s \ (69a) 

and 

a/ 2 =^/a VV 2 +^VV 3 . (72a) 

Here the constants A , A ' f M, A u A\ , A 2 and A 2 are functions of k 2 , 
k 3 , k 2 ', and k 3 ' and may be determined from eq 49, 51, 54, 70, 71, 73, 
and 74. For k 2} & 3 , k 2 ', and k z ' all equal to zero, we obtain 

A*=AJ =0.50000, 
Ai=Ai'=0.15000, 
A*=A/= 0.023936, 
M=L 

Now, the determination of the unknown quantities a , a u tj s and 
&o', <h') I* from the above six relations, 47a to 72a, is an easy task 
once a specific value is assigned to r. Separate determinations were 
made by selecting nine values for r, namely, r 2 =0.00, 0.01, 0.10, 0.3162, 
1, 3.162, 10.0, 100, and oo. The results are shown in table 5. 
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Table 5. — First approximate values of the constants ao, a u a\ ', 77, and 17/ 

7- 2 =(mV)/(mp) 



r 2 


a 1 


ai', ao 


or ao' 


V* 


«•' 


H. 


w 


0.00 
.01 




0. 6038 
.5402 


1.000 
0. 9285 




3.309 
3.433 




0. 995 
.959 


0. 0540 


2.611 


2.552 


.10 


.1388 


.4390 


.8085 


2.744 


3.678 


2.582 


.893 


.3162 


.2031 


.3614 


.7100 


2.845 


3.923 


2.595 


.838 


1.0 


.2727 


.2724 


.5910 


2.995 


4.331 


2.616 


.772 


3.162 


.3397 


.1911 


.4646 


3.146 


4.851 


2.630 


.679 


10 


.3916 


.1238 


.3479 


3.324 


5.605 


2.672 


.590 


100 


.4494 


.0449 


.1773 


3.655 


7.844 


2.750 


.423 


00 


.4839 




.000 


4.124 




2.884 





The expressions for H and H' from eq 43 and 47, respectively, 
wherein the constants k and ¥ are put equal to zero and the values 
of the basic constants are taken from table 5, constitute the first 
approximate solution of the laminar boundary layers at the interface 
of two liquids, one of which is still and the other moving. More 
specifically, in the upper liquid 



and in the lower liquid 



dH 

dr) 



=a +a 1 rj s G l (rj/r] s ) ) 



JTJt 

-j-r=a ' + a 1 'r) s 'G l (r J 'l7 1 /). 



The numerical values of F , F h and 6?i are given in tables 3 and 4. 

It may be of interest to note that the case r 2 = 00 can be physically 
interpreted as representing the case in which the lower liquid has 
an infinite viscosity, which is obviously the case for a liquid flowing 
over a plate. The solution obtained here for r 2 = 00 is identical with 
Pohlhausen's fourth approximation. In the present analysis this con- 
stitutes only the first approximate solution. 

VII. SECOND APPROXIMATE SOLUTION OF THE LAMINAR 
BOUNDARY LAYERS 

Inserting the first approximate values of a , a 1} and a ', d\ from 
table 5 in eq 78, 79, 80, and 81, we obtain a set of values for k 2 , k Z) 
k 2 ', and k z ', after discarding the squares and the higher powers of 
these constants. The results are given in table 6 for the various 
values of r 2 . 

Table 6. — Second approximate values of the expansion constants hi, £3, k 2 ' , and k$ 

r2=(Vp')/(MP) 



ri 


fej 


*a 


w 


h' 


M 


0.00 






-0. 2600 


+0. 0382 




.01 


-0. 1354 


-0. 0278 


-. 2600 


+. 0382 


1. 0681 


.10 


-.1259 


-.0328 


-.2600 


+. 0382 


1. 0737 


.3162 


-. 1162 


-.0380 


-. 2600 


-h 0382 


1. 0816 


1.0 


-.1037 


-.0446 


-. 2600 


+. 0382 


1. 0869 


3.162 


-.0897 


-.0522 


-. 2600 


+. 0382 


1. 0956 


10 


-.0734 


-. 0608 


-. 2600 


+. 0382 


1. 1057 


100 


-.0429 


-.0771 


-.2600 


+. 0382 


1. 1252 


00 


.0000 


-. 1000 
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With these, corresponding values of the coefficients Ap, A', A u 
Ai', A 2 , A 2 ' , and M are computed. The results are shown in table 7. 
With the values of the coefficients thus known, the second approxi- 
mate values of the basic constants are determined by the same 
equations as for the first approximate solution, that is, eq 48a to 72a, 
inclusive. The results of this determination are given in table 8. 

Expressions for H and H' from eq 43 and 47, respectively, where 
the expansion constants are given in table 6 and the basic constants 
are given in table 8, constitute the second approximate solution of 
the laminar layers boundary problem for the interface between two 
liquids. 

Table 7.- — Second approximate values of the auxiliary constants A , A h A 2 , ^V, ^i', 

and A%* 

r2=( M V)/(Mp) 



r 2 


.lo 


At 


A% 


Ao' 


Ai 


Ai 


0. 00 
.01 








0. 3806 
.3806 


0. 08562 
. 08562 


0.01336 
. 01336 


0. 4139 


0. 1110 


0. 01830 


.10 


.4164 


.1131 


. 01871 


.3806 


. 08562 


. 01336 


.3162 


.4190 


.1155 


.01917 


.3806 


. 08562 


. 01336 


1.00 


.4224 


.1181 


. 01970 


.3806 


. 08562 


.01336 


3. 162 


.4260- 


.1212 


. 02038 


.3806 


. 08562 


.01336 


10 


.4304 


. 1250 


. 02107 


.3806 


. 08562 


. 01336 


100 


.4382 
.4600 


.1321 
.1427 


. 02245 
. 02460 


.3806 


. 08562 


. 01336 



Table 8.- — Second approximate values of the basic co7istants a , «i, «i', is, and r) Q 

r2=( Al 'p')/(/xp) 



r 3 


a\ 


ai 


ao or ao' 


Vs 


ni 


H» 


Hi 


0.00 

.01 




-0. 5906 
-.5329 


1. 0000 
0. 9282 




4.449 
4.576 


3. 012 
3.004 


1.041 
0. 993 


0. 0567 


3.069 


.10 


.1456 


-. 4287 


.8082 


3.163 


4.953 


2.991 


.936 


.3162 


.2144 


-.3527 


.7088 


3.241 


5.280 


2.990 


.875 


1.0 


.2899 


-.2668 


.5888 


3.358 


5.793 


2.991 


.798 


3.162 


.3611 


-.1854 


.4619 


3.498 


6.545 


2.996 


.707 


10 


.4174 


-.1194 


.3446 


3.648 


7.582 


3.016 


.602 


100 


.4817 


-. 0428 


.1739 


3.913 


10. 672 


3.055 


.434 


oo 


.5196 




.0000 


4.277 









VIII. QUALITY OF THE APPROXIMATIONS 

The quality of the approximations may be studied conveniently by 
substituting the approximate expressions for H and H f and for their 
derivatives into the differential equations, eq 12 and 15, and then 
considering the magnitude of the remainders. Let the remainders 
be denoted by Ai and Ai' : 



AZ= 



and 

If the ratios 






Ai/((fH/rf7j 3 ) mttX and AL'Kd'H'/dr,' 
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become smaller in the successive approximations, the magnitude of 
decrease is a measure of the improvements in the approximations. 
We have computed these remainders only for the two extreme 

0.24 



0.18 






«r° o.o6 



0.00 



-0.06 
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Figure 2 — Remainders of the first and of the second approximations for the case 
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Figure 3 — Remainders of the first and of the second approximations for the case 
r*=0. 

cases represented by r 2 =co and r 2 —0. The values of the remainders 
together with the values of the third derivatives of H and H', are 
given in figures 2 and 3, respectively. The case of a liquid of finite 



- — _ 
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viscosity moving over a liquid of infinite viscosity is represented by 
r 2 =oo. As mentioned before, this is mathematically identical with 
the Blasius problem of flow over a plate. In examining the curves 
in figure 2, we see that the amplitude of the oscillation of the remainder 
in the second approximate solution is only about half as large as in the 
first approximation, indicating that the improvement in the second 
approximate solution is appreciable. In the second approximate 
solution, the maximum value of the remainder is 0.027; the maximum 
value of the third derivative, 0.25. It is thus seen that the remainders 
are relatively small. In this same figure we give also the curve of the 
third derivative of H from the exact solution by Dryden [5]. The 
agreement between the second approximate solution and the exact 
solution is rather close. The velocity gradient at y=0, from the 
second approximate solution, is 

\j ^=0.3307 ■mW. 

In the exact solution the corresponding coefficient has the value 
0.3320, thus indicating that the error in the value of this quantity 
from the second approximate solution is only 0.39 percent. 

A liquid of infinite viscosity moving over a liquid of finite viscosity 
is represented by r 2 =0. Physically, it would seem that this is the 
case of a solid surface moving in a still liquid. The curves in figure 3 
also show that in this case the remainder of the second approxi- 
mate solution has been reduced and has become about half as large as 
the corresponding quantity in the first approximate solution. But the 
accuracy obtained is not as satisfactory as in the case r 2 = oo ? since 
the remainders are relatively larger. Apparently to obtain the same 
degree of accuracy as was obtained in the case r 2 = co } and with the 
second approximation it will be necessary to introduce the expansion 
functions kJVi and k 5 N 5 . 

IX. PHYSICALLY SIGNIFICANT QUANTITIES OF THE 

LAYERS 

Among the quantities of the laminar layers at the interface of the 
two liquids, one at rest and the other moving over it, the following 
are of practical significance: (1) The intensity of the viscous stresses 
at the surface of separation, (2) the thickness of the layers, and (3) 
the distribution of the velocities in the layers. The computation of 
these quantities will be made, using the second approximate solution 
of the analysis. 

1. LOCAL STRESS AT THE SURFACE OF SEPARATION 

The stress in question is determined from the velocity gradient of 
the upper liquid at the surface of separation: 

From eq 19 

U/d 2 H\ 



. ... — 
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or 



T -§rA^ 2 /o2 U2 ' 



Substituting from eq 14, 

*-<f).(f r^- 

If we use the velocity gradient of the lower liquid, we obtain in the 
same manner, 



Putting 
and 



F JcPH'\(Ux\-* P ' U2 
-v , 2(^) o =V2a 1 (l+/c 3 )= > s J 



we have the local stresses 



or 






(82) 



(83) 



The computed values of s and s' as functions of r 2 are given in 
table 9. 

Table 9. — Values of the constants s and s' in the expressions of the local stress 
at the interface, eq 82 and 88 

r 2 =(MV)/(/«p) 



r 2 


S 


s' 


0.00 




0. 8672 


.01 


0. 0782 


.7825 


.10 


.1992 


.6295 


. 3162 


.2917 


. 5178 


1.0 


. 3917 


.3918 


3. 162 


.4840 


.2722 


10 


. 5543 


.1753 


100 


.6287 


.0628 


CO 


.6612 







2. THICKNESSES OF THE LAYERS 

Since the velocities in the laminar boundary layers approach the 
limiting values in an asymptotic manner, the definition of the thick- 
ness of the layers becomes a matter of convention. Of particular 
utility are the following definitions. See figure 4. 



>r 



2J o (U-u)dy=(U-u )5 1> 



(84) 
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Figure 4. — The gradient thicknesses of the interfacial laminar boundary layers. 
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and 

2 f * u'dy^iiffit', (85) 

the first applying to the upper moving liquid and the second to the 
lower still liquid. It will be appropriate to speak of 5i and of 5/ as 
the gradient thicknesses. Like the displacement thicknesses, the 
values of 6 X and 6/ would not be affected very much by the approxi- 
mations, and this obviously, can not be said of y s and y/ . Using the 
relations 

y=rfi, y'=V«', 

8=- y f2vx/U i 

u=UdH/drj, 

we find from eq 84 and 85 that 



S' = ij2v'x/U, 
u' = UdH'/d v ', 



vx 



and 



Putting 



v=^fVrr 



n^2^2 r f. 



& n'=2J2^, 
-do a 



we have for the gradient thicknesses, 



$i=nJ 



and 5/= ft' 






(86) 



The values of n and n' are given in table 10 as functions of r 2 . We 
obtain for the boundary thicknesses, introducing the expressions 



Vvx , , \v'x 



(87) 



The values of the coefficients m and m' as functions of r 2 are given in 
table 10. 

Table 10. — Values of the constants n and n' in the expressions of the gradient thick- 
ness and m and m' in the expressions of the boundary thicknesses of the two lam- 
inar layers t eq 86 and 87 

r*-<nV)/(«0 



ra 


n 


n' 


m 


m' 


0.00 




2.943 




6.291 


.01 


2.270 


3.026 


4. 339 


6.470 


.10 


2.344 


3.304 


4.472 


7.003 


.3162 


2.428 


3.493 


4.582 


7.466 


1.0 


2.531 


3.831 


4.748 


8.191 


3.162 


2.664 


4.328 


4.946 


9.255 


10 


2.813 


4.940 


5.158 


10. 721 


100 


3.071 


7.063 


5.533 


15. 090 


00 


3.455 




6.047 
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3. DISTRIBUTION OF THE VELOCITIES IN THE LAYERS 

Those velocities are determined from eq 42 and 46. The results of 
the computations for the upper liquid are given in table 11, where 
u/D is tabulated as a function of y/y s for various values of r 2 . Similar 
results for the lower liquid are given in table 12. 

Table 11. — Velocity distribution in the inter jacial laminar layer of the upper, 

moving liquid 

The boundary thickness y„ is fdven by eq 87 

r 2 =OV)/(/*p) 



X 


















\ r2 

V/Vb \ 


0.010 


0.100 


0.316 


1.000 


3.16 


10 


100 


CO 


0.0 


v/U 


0. 9282 


0. 8082 


0. 7088 


0. 5888 


0. 4619 


0. 3446 


0. 1739 


0.0000 


.1 


.9449 


. 8521 


.7747 


.6807 


.5808 


. 4862 


. 3466 


.1994 


.2 


.9602 


.8929 


.8363 


.7667 


.6917 


.6204 


.5122 


.3935 


.3 


. 9735 


.9279 


.8892 


. 8412 


.7891 


.7387 


.6609 


.5723 


.4 


.9838 


.9557 


.9316 


. 9015 


.8685 


.8359 


.7850 


. 7271 


.5 


.9913 


.9758 


.9626 


.9457 


.9273 


.9089 


.8797 


.8440 


.6 


.9961 


.9889 


.9826 


.9747 


.9660 


. 9573 


.9432 


.9252 


.7 


.9988 


.9962 


.9933 


.9910 


.9878 


.9848 


.9801 


.9727 


.8 


.9998 


.9991 


.9985 


.9980 


.9973 


.9966 


.9960 


.9938 


.9 


.9999 


.9999 


.9998 


.9998 


.9997 


.9996 


.9996 


.9995 


1.0 


1.0000 


1. 0000 


1. 0000 


1. 0000 


1.0000 


1. 0000 


1. 0000 


1. 0000 



Table 12. — Velocity distribution in the inter jacial laminar layer of the lower, 

still liquid 

The boundary thickness y t ' is given by oq 87 

r 2 =0xV)/(w>) 



N 




















0.00 


0.010 


0.100 


0.316 


1.000 


3.16 


10 


100 


0.0 








u' 


tu 








1. 0000 


0.9282 


0. 8082 


0. 7088 


0. 5888 


0. 4619 


0. 3446 


0. 1739 


.1 


.7337 


.6811 


.5930 


.5200 


.4320 


.3389 


.2527 


.1276 


.2 


.4980 


.4624 


.4023 


.3530 


.2932 


.2301 


.1716 


.0866 


.3 


.3121 


.2897 


.2524 


.2213 


.1838 


.1441 


.1075 


.0542 


.4 


.1765 


.1639 


.1427 


.1250 


.1039 


.0814 


. 0608 


.0306 


.5 


.0886 


.0822 


.0717 


.0628 


. 0522 


.0404 


.0305 


.0154 


.6 


.0369 


. 0343 


.0298 


.0262 


.0217 


.0171 


.0131 


.0064 


.7 


.0119 


.0112 


.0098 


.0092 


.0070 


.0055 


.0041 


.0020 


.8 


.0026 


.0025 


.0022 


.0021 


.0016 


.0011 


.0009 


.0004 


.9 


.0007 


.0006 


.0005 


.0004 


.0003 


.0002 


.0001 


.0001 


1.0 


.0000 


.0000 


.0000 


.0000 


.0000 


.0000 


.0000 


.0000 



A more suggestive picture of the velocity distribution is obtained 
by plotting uj U against y/y s , y now being the distance from the surface 
of separation to a point, above or below the surface, and y s being the 
boundar}^ thickness of the layer situated in the upper liquid. Figure 
5 gives the velocities for identical liquids, p— j»f and p=p'. Figure 6 
gives the velocities for v' = 7.85v and p / = 1.14p. This case will be 
realized if we take water for the upper moving liquid and aqueous 
solution of 56 percent glycerol for the lower, still liquid, both being at 
20° C. These two figures demonstrate the very large influence that the 
increased viscosity of the lower liquid has on the velocity distribution. 
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Figure 5. — Interfacial laminar velocity distribution for the case r 2 —l, 
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Figure 6. — Interfacial laminar velocity distribution for the case r 2 — 10. 

X. TRANSVERSAL VELOCITIES OF THE LAYERS 

The expressions of the transversal velocities, v and v' , in terms of 
H and v, may be derived conveniently from the condition of continuity 
eq 4. Since at the interface, v and v' vanish, we have 
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This yields, using eq 13, 14, and 18, the relation 



j C"cPH , 



Integrating by parts and noting that H=0 when 77=0, we obtain 

?-f«-* « 

which is the expression for the transversal velocities at the point 

(1?, *). 

The corresponding expression for the lower fluid is 

b'v'_dH' 

-V~W V ~ a • (89) 

Thus, at the points of the lower liquid and in the region outside of the 
layer, 7i'^>rj s ', there is a weak current moving normally upward 
toward the interface. The strength of the current is given by 

In a sense, the boundary layer of the upper liquid acts as a pump, 
raising the small portions of the lower liquid to the level of the inter- 
face, then causing these portions to move horizontally. So long as 
these motions as required by the theory are not interfered with, the 
steadiness of the layer will be assured. 
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